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Background: At energies near the Coulomb barrier, capture reactions in heavy-ion collisions result either in 
fusion or in quasifission. The former produces a compound nucleus in statistical equilibrium, while the second 
leads to a reseparation of the fragments after partial mass equilibration without formation of a compound nucleus. 
Extracting the compound nucleus formation probability is crucial to predict superheavy-element formation cross- 
sections. It requires a good knowledge of the fragment angular distribution which itself depends on quantities 
such as moments of inertia and excitation energies which have so far been somewhat arbitrary for the quasifission 
contribution. 

Methods: We investigate the evolution of the nuclear density in time-dependent Hartree-Fock (TDHF) calcu¬ 
lations leading to quasifission. Our main goal is to extract ingredients of the formula used in the analysis of 
experimental angular distributions. These include the moment-of-inertia and temperature. We study the depen¬ 
dence of these quantities on various initial conditions of the reaction process. 

Results: The evolution of the moment of inertia is clearly non-trivial and depends strongly on the characteristics 
of the collision. The temperature rises quickly when the kinetic energy is transformed into internal excitation. 
Then, it rises slowly during mass transfer. 

Conclusions: Fully microscopic theories are useful to predict the complex evolution of quantities required in 
macroscopic models of quasifission. 

PACS numbers: 21 . 60 .-n, 21 . 60 .Jz 


I. INTRODUCTION 

The creation of new elements is one of the most novel 
and challenging research areas of nuclear physics [1-4]. 
The search for a region of the nuclear chart that can 
sustain the so called superheavy elements (SHE) has led 
to intense experimental activity resulting in the discov¬ 
ery and confirmation of elements with atomic numbers 
as large as Z = 117 [5-7]. The theoretically predicted 
island of stability in the SHE region of the nuclear chart 
is the result of new proton and neutron shell-closures, 
whose location is not precisely known [8-10]. The exper¬ 
iments to discover these new elements are notoriously dif¬ 
ficult, with fusion evaporation residue (ER) cross-section 
in pico-barns. This cross-section is commonly expressed 
in the product form 

J max 

^ER — ^ ^ ^cap {Ec.m.,J)PcN{E^, J)Wsur{E^, J) ( 1 ) 

J=0 

where (Jcap{Ec.m.^ J) is the capture cross section at center 
of mass energy Ec.m. and spin J. Fcn is the probability 
that the composite system fuses into a compound nucleus 
(CN) rather than breaks up via quasifission. Wsur is the 
survival probability of the fused system against fission. 
Eor light and medium mass systems the capture cross- 
section may be considered to be the same as that for 
complete fusion. However, for heavy systems leading to 
superheavy formation, the ER cross-section is dramati¬ 
cally reduced due to the quasi-fission and fusion-fission 
processes [11, 12], thus making the capture cross-section 
to be essentially the sum of these two cross-sections. 


The fusion process implies a transition from a dinuclear 
configuration - accompanied by particle exchange during 
the dynamical process - to a single-center compound-like 
configuration of the composite system. Most dynamical 
models [13-21] argue that, for heavy systems, a dinuclear 
complex is formed initially and the barrier structure and 
the excitation energy of this precompound system deter¬ 
mine its survival to breaking up via quasi-fission. Eur- 
thermore, if the nucleus survives this initial state and 
evolves to a compound system, it can still fission due to 
its excitation. 

Among the three stages towards ER given in Eq. (1), 
the determination of Fcn contains the most uncer¬ 
tainty which can be as much as 1—2 orders of magni¬ 
tude [22, 23]. Experimentally, Pew can be extracted 
from the measurement of fusion-evaporation residue 
cross-sections [24, 25]. However, these cross-sections be¬ 
come very small for heavy systems and the extraction of 
PcN is uncertain for such systems. Information on Pcn 
can also be obtained by comparing the width of fragment 
mass distribution with the width expected in the case of 
pure fusion-fission [26, 27]. This approach, however, only 
provides an upper limit for Pcn as it assumes that only 
fusion-fission produces symmetric fragments [26]. 

An alternative approach involves the analysis of the 
fragment angular distribution. Eor instance, assuming 
that the quasifission process, as fusion-fission, is sta¬ 
tistical, the critical angular momentum Jcn between 
fusion-fission and quasifission can be adjusted to repro¬ 
duce the experimental angular distribution of symmetric 
fragments [23, 28, 29]. Although this approach has been 
widely used, it could have limitations due to the fact 
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that quasifission in not a statistical decay but a dynam¬ 
ical process, and that fusion-fission could be asymmet¬ 
ric due to, e.g., late-chance fusion-fission at low excita¬ 
tion energies [30] as well as shell structure of pre-scission 
configurations [31]. Correlations between mass and scat¬ 
tering angle of the fragments have also been measured 
extensively [27, 32-35]. In particular, they can be used 
to disentangle fast quasi-fission processes (few zeptosec- 
onds) to longer reaction mechanisms associated with con¬ 
tact times between the fragments exceeding 10 — 20 zs 
(i.e., long-time quasifission and fusion-fission). However, 
the analysis of such fragment mass-angle distributions 
[32, 33], as well as statistical descriptions of fragment an¬ 
gular distributions [28] , require external parameters such 
as moment of inertia and temperature which are chosen 
somewhat arbitrarily [23, 28, 36]. It is therefore impor¬ 
tant to provide realistic evaluation of these parameters, in 
particular for the quasifission mechanism which requires 
a description of the complex nuclear dynamics. 

Dynamical microscopic approaches are a standard tool 
to extract macroscopic properties in heavy-ion collisions 
[37-41]. In particular, the time-dependent Hartree-Fock 
theory [42], has been recognized for its realistic descrip¬ 
tion of several low-energy nuclear reaction mechanisms 
[43, 44]. It has been recently utilized for studying the 
dynamics of quasifission [27, 35, 45, 46] and scission [47- 
49] . The study of quasifission is showing a great promise 
to provide insight based on very favorable comparisons 
with experimental data [27, 35]. As a dynamical micro¬ 
scopic theory, TDHF provides us with the complete shape 
evolution of the nuclear densities which can be used to 
compute the time evolution of deformation and inertia 
parameters. Using an extension to TDHF, the so-called 
density constraint TDHF (DC-TDHF) approach [37], it 
is also possible to compute the excitation energy in the 
fragments. 

In this manuscript we focus our discussion on the ex¬ 
traction of the time evolution of the moment of inertia 
and of the excitation energy (temperature). These are 
indeed the relevant quantities for fragment mass-angle 
distributions and angular distribution analyses. These 
quantities are computed for collisions between calcium 
isotopes and actinides for which the TDHF approach has 
been shown to provide a deep insight into the quasifission 
reaction mechanisms [35, 45, 46]. 


II. FORMALISM 

A. TDHF and DC-TDHF approaches 

The theoretical formalism for the microscopic descrip¬ 
tion of complex many-body quantum systems and the un¬ 
derstanding of the nuclear interactions that result in self¬ 
bound, composite nuclei possessing the observed proper¬ 
ties are the underlying challenges for studying low energy 
nuclear physics. The Hartree-Fock approximation and 
its time-dependent generalization, the time-dependent 


Hartree-Fock theory, have provided a possible means to 
study the diverse phenomena observed in low energy nu¬ 
clear physics [43, 44]. In general modern TDHF calcula¬ 
tions provide a useful foundation for a fully microscopic 
many-body description of large amplitude collective mo¬ 
tion including collective surface vibrations and giant res¬ 
onances [50-63] nuclear reactions in the vicinity of the 
Coulomb barrier, such as fusion [38, 41, 64-71], deep- 
inelastic reactions and transfer [21, 72-77], and dynamics 
of (quasi)fission fragments [35, 45-49]. 

Despite its successes, the TDHF approach has impor¬ 
tant limitations. In particular, it assumes that the many- 
body state remains a single Slater determinant at all 
time. It describes the time-evolution of an independent 
particle system in a single mean-field corresponding to 
the dominant reaction channel. As a result, it induces 
a classical behavior of many-body observables. A well 
known example is that TDHF does not include tunneling 
of the many-body wave-function and, thus, it is unable to 
describe sub-barrier fusion. It also underestimates width 
of fragment mass and charge distributions in strongly dis¬ 
sipative collisions [72, 74, 78]. Thus, to obtain multiple 
reaction channels or widths of one-body observables one 
must in principle go beyond TDHF [74, 79-82]. 

In fact, different reaction outcomes can also be ob¬ 
tained from mean-field calculations if the initial state 
(defined by the center of mass energy Ec.m. and by the 
orbital angular momentum L of the collision) is best de¬ 
scribed by a superposition of independent (quasi)particle 
states and that each of these states can be assumed to 
evolve in its own mean-field. This is the case, for in¬ 
stance, if a collision partner is deformed in its intrinsic 
frame. In this case, each orientation of the deformed 
nucleus may encounter a different mean-field evolution 
[83, 84]. Such orientation dependence of reaction mech¬ 
anisms has been experimentally studied in quasifission 
with actinide targets [85-91] and confirmed in TDHF 
studies [27, 35, 44]. In particular, these studies have 
shown that collisions of a spherical projectile with the tip 
of prolately deformed actinides lead to fast quasifission 
(with contact time smaller than 10 zs) while collisions 
with the side of the actinide may induce longer contact 
times, larger mass transfer, and possible fusion. It is 
therefore common practice to investigate a subset of spe¬ 
cific orientations depending on the reaction mechanism 
one is interested to investigate. For instance, it is suffi¬ 
cient to study ’’side collisions” in order to investigate the 
competition between quasifission and fusion in collisions 
with actinide nuclei [45]. 

In recent years has it become numerically feasible to 
perform TDHF calculations on a 3D Cartesian grid with¬ 
out any symmetry restrictions and with much more accu¬ 
rate numerical methods [44, 54, 75, 92-94]. In addition, 
the quality of effective interactions has been substantially 
improved [95-98] . In order to overcome the lack of quan¬ 
tum tunneling preventing direct studies of sub-barrier fu¬ 
sion, the DC-TDHF method was developed to compute 
heavy-ion potentials [37] and excitation energies [99] di- 
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rectly from TDHF time-evolution. For instance, this 
method was applied to calculate capture cross-sections 
for hot and cold fusion reactions leading to superheavy 
element Z = 112 [67]. 


B. Fragment angular distributions 


component expression for fusion-fission and quasifission 
parts [28, 100-102], 


w{e) = ^Ff^\e)+ £ ( 2 ) 

t7=0 J=JcN 


Experimental analysis of the fragment angular distri¬ 
butions W{0) is commonly expressed in terms of a two- where 


_ (2J -1- l)^exp[-{J -1- l/2)^sin^0/4K^{a)]Jo[i{J + l/2)^sin^0/4K^{a)] 
^ ~ er/[(J+l/2)/(2X2(a))V2] 


and a = FF (fusion-fission) or QF (quasi-fission). Here, 
JcN defines the boundary between fusion-fission and 
quasifission, assuming a sharp cutoff between the angu¬ 
lar momentum distributions of each mechanism. The de¬ 
tailed definition of various mathematical functions can 
be found in Refs. [23, 29]. 

The quantum number K is known to play an important 
role in fission [103]. The latter is defined as the projection 
of the total angular momentum along the deformation 
axis. In the Transition State Model (TSM) [103], the 
characteristics of the fission fragments are determined by 
the K distribution at scission. The argument Kq entering 
Eq. 3 is the width of this distribution which is assumed 
to be Gaussian. It obeys 

= T%ff/h^ , (4) 

where the effective moment of inertia, /, is computed 
from the moments of inertia for rotations around the axis 
parallel and perpendicular to the principal deformation 
axis 


and T is the nuclear temperature at the saddle point. 
The physical parameters of the fusion-fission part are rel¬ 
atively well known from the liquid-drop model [104, 105]. 
In contrast, the quasifission process never reaches sta¬ 
tistical equilibrium. In principle, it has to be treated 
dynamically, while equation 3 is based on a statistical 
approximation. In addition, the usual choice for the nu¬ 
clear moment of inertia for the quasifission component, 
^o/^e// = 1-5 [23, 28, 36], is somewhat arbitrary. Here, 
7^0 is the moment of inertia of an equivalent spherical 
nucleus. 

In order to illustrate the sensitivity of the fragment an¬ 
gular distribution with the chosen value of the moment of 
inertia, we have plotted the ratio 1F(6>)/1F(90°) in Eig- 
ure 1 with T^o/^e// = 1-1 and 1.5. Here, we used the 
value of JcAT = 46 obtained directly from TDHE cal¬ 
culations of quasifission for this system. Similarly, tem¬ 
perature was taken to be T = 1.2 MeV. The deviation 



FIG. 1. (Color online) Effect of the moment of inertia on the 
angular distribution of the fragments computed with Eqs. 2 
and 3. Both curves were calculated using Jcn = 46 and 
1.2 MeV. 


is mostly visible at most forward and backward angles, 
where it can reach up to 20%. It is interesting to note 
that if we assume that the red curve of Eig. 1 is the actual 
measured quantity and try to fit it by varying the Jcn 
value of the black curve we obtain a value oi Jcn = 37. 
This would then be the error for using the wrong Jcn 
value. 

In the following section we outline the extraction of 
these ingredients directly from TDHE time-evolution of 
collisions resulting in quasifission. 


III. RESULTS 

The feasibility of using TDHE for quasifission has only 
been recognized recently [35, 44, 45]. By virtue of long 
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contact-times for quasifission and the energy, orientation 
and impact parameter dependencies, these calculations 
require extremely long CPU times and numerical accu¬ 
racy [92-94, 106]. During the collision process, the nu¬ 
clear densities, as described by TDHF time-evolution, un¬ 
dergo complicated shape changes including vibration and 
rotation. Such evolutions finally lead to two separated 
final fragments identified as quasifission due to the long 
contact-time for the reaction as well as the mass/charge 
of the fragments [35, 45]. 



FIG. 2. (Color online) Quasi-fission in the reaction 
^^Ca+^^^Bk at Fc.m. = 218 MeV with impact parameter 
b = 2.0 fm. Shown is a contour plot of the time evolution 
of the mass density. The actual numerical box is larger than 
the one shown in the frames. 

In Fig. 2 we show a few time snapshots of the evolv¬ 
ing mass density for the QF reaction of the ^^Ca+^^^Bk 
system at F^c.m. = 218 MeV with impact parameter 
6 = 2.0 fm. The times are given in zeptoseconds {1 zs = 
10“^^s). The very large elongation of the separating 
fragments is noteworthy. The initial orientation of the 
deformed ^^^Bk was chosen such that the spherical ^^Ca 
nucleus collides with the side of the actinide nucleus. In 
the following we will refer to this as the “side” orienta¬ 
tion. 


A. Incoming and outgoing potentials 

As quasifission is part of the capture process, it is 
therefore possible to compute the nucleus-nucleus poten¬ 
tial down to the inside barrier region with the DC-TDHF 
technique. The latter is represented in Fig. 3 for the 
incoming channel of a "^^Ca+^^^U central collision. Af¬ 
ter contact, the two fragments encounter a significant 
mass transfer toward symmetry. The outgoing potential 
is therefore different than the incoming one. In fact, the 
outgoing potential does not exhibit any barrier, leading 
to a re-separation of the fragments. 



FIG. 3. (Golor online) DG-TDHF nucleus-nucleus potential 
for the ^°Ga+^^®U central collision at Ec.m. =211 MeV and 
with the side orientation. Both the potentials for the incoming 
and outgoing channels are shown. 


B. Moment of inertia 

The first collective observable of interest for fission and 
quasi-fission (both dynamical and statistical) studies is 
the moment of inertia of the system. The proper way to 
calculate the moment-of-inertia for such time-dependent 
densities (particularly for non-zero impact parameters) 
is to directly diagonalize the moment-of-inertia tensor 
represented by a 3 x 3 matrix with elements 

%j{t)/m = J d^r - XiXj) , (6) 

where p is the local number-density calculated from 
TDHF evolution, m is the nucleon mass, and ^^= 1 ^ 2,3 de¬ 
note the Cartesian coordinates. The TDHF calculations 
are done in three-dimensional Cartesian geometry [93]. 
Numerical diagonalization the matrix ^ gives three eigen¬ 
values. One eigenvalue corresponds to the moment-of- 
inertia ^11 for the nuclear system rotating about the prin¬ 
cipal axis. The other two eigenvalues define the moments 
of inertia for rotations about axes perpendicular to the 
principal axis. Naturally, for triaxial density distribu¬ 
tions, the two perpendicular components are not exactly 
equal but for practical calculations they are close enough 
and always larger than the parallel component. We thus 
use a single average value for these moments of inertia 
denoted by ^_l- 

Using the time-dependent moment-of-inertia obtained 
from the TDHF collision one can calculate the so-called 
effective moment-of-inertia defined in Eq. (5). It is stan¬ 
dard to compute the effective moment of inertia rela¬ 
tive to a spherical system using the mass independent 
quantity where ^0 is the moment-of-inertia 
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FIG. 4. TDHF results showing the time-dependence of the ra¬ 
tio ^o/^eff for the system at Ec.m. = 218 MeV 

and for impact parameters ranging from 6 = 0.5 to 3 fm. The 
arrows indicate the scission point of each trajectory. 


FIG. 5. TDHF results showing the time-dependence of the 
ratio ^o/^eff for the ^°Ga + system at energies Ec.m. = 
208 and 211 MeV for zero impact parameter. 


of a spherical nucleus with the same number of nucle¬ 
ons [23, 29]. The expression for the moment-of-inertia 
for a rigid sphere is given by ^o/m = 2ARqI5^ where Rq 
can be chosen as Rq = 1.22bA^^^ fm [29]. 


1. Role of impact parameter 

The moment-of-inertia ratio calculated for the 
^^Ca + non-central collisions at Ec.m. = 218 MeV 

is shown in Fig. 4. These trajectories all lead to quasi¬ 
fission with a large mass transfer and orbiting before the 
separation of the fragments. This is illustrated in Fig. 2 
for the impact parameter 6 = 2 fm case. For this system 
using the value of h‘^/m = 41.471 MeV-fm^ and total 
mass number A = 297, we get % = 191.359 h‘^'MeV~^. 
At the point of final touching configuration the moment- 
of-inertia ratios are in the range 1.4 — 1.8, suggesting a 
relatively strong impact parameter dependence. Impact 
parameters smaller than 6 = 0.5 fm lead to contact time 
between the fragments exceeding 35 zs. We can therefore 
consider that the nuclei have fused and would ultimately 
form a compound nucleus. 


2. Role of center of mass energy 

Next we consider central collisions of the ^^Ca+^^^U 
systems. Similar to the ^^^Bk case, exhibits a strong 
prolate deformation and its alignment with respect to the 
collision axis changes the quasifission characteristics [35] . 
The value of is 171.394 h^'MeV~^ for A = 278. 

Figure 5 shows the time-evolution of the moment-of- 


inertia ratio at two different center of mass energies for 
central collisions with the side orientation. The two cal¬ 
culations have an energy difference of only 3 MeV in order 
to compare reactions which lead to similar mass asym¬ 
metry of the fragments. It is interesting to observe that 
up to ^ 8 zs, the two evolutions of the moment of inertia 
are almost identical. This is likely to be due to the fact 
that the translational kinetic energy is rapidly dissipated 
in both cases, leading to a system with similar compact¬ 
ness but different excitation energies. For longer times, 
we observe a faster re-separation of the fragments at the 
highest center of mass energy. However, the final value 
of the ratio — 1-2 is approximately the same 

for both energies. Note that this value is again differ¬ 
ent from the assumed ratio of 1.5, but in this case it is 
significantly smaller than in the ^^Ca + ^^^Bk system, 
indicating that the ^^Ca + system leads to more 

compact configurations. 


3. Role of the orientation 

Figure 6 shows the dependence of the moment-of- 
inertia ratio on the orientation of the ^38^ nucleus for 
the "^^Ca+^^^U system (^o = 179.693 /i^-MeV“^ for 
A = 286). The shorter contact time for the tip orien¬ 
tation (elongation axis of parallel to the collision 

axis) is evident from the figure. We also observe that at 
the initial contact and for most of the evolution, the ra¬ 
tios are significantly different, originating from the very 
different initial density configurations. However, as the 
two reactions approach the scission point the ratios be¬ 
come relatively close to one another: The neck breaks at 
^o/^eff ^ 1.1 — 1.2. These values are, again, signifi- 
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cantly smaller than the assumed value of 1.5. 



FIG. 6. TDHF results showing the time-dependence of the 
ratio ^o/^eff for the system at energy Ec.m. = 

203 MeV for zero impact parameter and two orientations of 
the nucleus. 


C. Collective dynamics and temperature 

In this Section we study the evolution of the TDHF 
density relative to the dynamical potential energy sur¬ 
face (PES) using the DC-TDHF method. In the DC- 
TDHF method a parallel static calculation is performed 
at given time intervals which constrains the instanta¬ 
neous TDHF density and finds the corresponding min¬ 
imum energy state. This allows us to trace the TDHF 
dynamical trajectory on the PES, albeit restricted to the 
shapes calculated in TDHE evolution. At the same time 
these calculations provide the dynamical excitation en¬ 
ergy, E'*(t) [99]. The DC-TDHE calculations for heavy- 
systems are extremely compute extensive and the calcu¬ 
lation of the trajectories performed here took about two 
months of computing on a 16 processor modern worksta¬ 
tion utilizing all processors. 

The temperature used in Eq. (4) is taken to be the 
temperature at the saddle point of a fissioning system [23, 
36] . Since this is adopted from modeling of fission it may 
not be appropriate for quasifission, which does not really 
have a saddle point in the sense of ordinary fission and 
it is a fully dynamical process. 

We have computed the dynamical temperature of the 
system during the quasifission path using the calculated 
excitation energy, as T{t) = i/E*(t)/(A/8.5) 

in MeV. Eigure 7(a) shows the TDHE results for the 
^^Ca + central collision at Ec.m. = 211 MeV and 
side orientation. We see that the temperature rapidly 
rises during the initial overlap phase of the collision and 
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FIG. 7. TDHF results showing the time-dependence of the 
temperature for the ^°Ga + system at energy Ec.m. = 

211 MeV for zero impact parameter and side orientation of 
the nucleus. 


reaches a value of about 1.0 MeV. After this the tempera¬ 
ture rises more slowly being in the range of 1.25 —1.5 MeV 
along the scission path. 

While the calculation of the full temperature dynam¬ 
ics is very computationally expensive one can calculate 
the temperature if the dynamical density is stored at a 
chosen time. This would then allow us to obtain all the 
ingredients of the parameter Kq of Eq. (4). One problem 
is that so far we do not have a convention for knowing 
exactly at what point one is supposed to evaluate the 
ingredients of Kq. The closest configuration to a fission 
saddle point would be the most compact shape formed 
just after full dissipation of the translational kinetic en¬ 
ergy. In this case, this would correspond to the end of 
the first, fast rise of T*(t), with a value of ^ 1 MeV. 


IV. SUMMARY 

The fully microscopic TDHE theory has shown itself 
to be rich in nuclear phenomena and continues to stim¬ 
ulate our understanding of nuclear dynamics. We have 
used the TDHE theory to study evolution of the nuclear 
density for reactions resulting in quasifission with a focus 
on the ingredients that are used to analyze experimen¬ 
tal angular distributions to calculate the fusion proba¬ 
bility PcN- We show that a number of useful quantities 
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can be obtained from TDHF dynamics rather than utiliz¬ 
ing models that may not be appropriate for quasifission. 
Among these are the moment-of-inertia parallel and per¬ 
pendicular to the symmetry axis, temperature through 
the calculation of the dynamical excitation energy, as well 
as other quantities like the rotational energy. In addition 
TDHF can tell us the dependence of these variables on 
impact parameter, energy, and structure. 
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